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Abstract 



We develop a new approach to Baxter Q-operators by relating them to the theory of Yangians, 
which are the simplest examples for quantum groups. Here we open up a new chapter in this 
theory and study certain degenerate solutions of the Yang-Baxter equation connected with har- 
monic oscillator algebras. These infinite-state solutions of the Yang-Baxter equation serve as 
elementary, "partonic" building blocks for other solutions via the standard fusion procedure. As 
a first example of the method we consider sl(n) compact spin chains and derive the full hierarchy 
of operatorial functional equations for all related commuting transfer matrices and Q-operators. 
This leads to a systematic and transparent solution of these chains, where the nested Bethe 
equations are derived in an entirely algebraic fashion, without any reference to the traditional 
Bethe ansatz techniques. 



1 Introduction and Overview 



The method of functional relations and commuting transfer matrices, introduced by Baxter in his 
seminal paper [1] on the exact solution of the eight-vertex model, plays a fundamental role in the 
theory of integrable quantum systems. It is based on an explicit algebraic construction of transfer 
matrices, which gives a priori knowledge about the analytic properties of their eigenvalues. 
A central part within this method involvess the so-called Q-operators. These operators are 
distinguished by the fact that zeroes of their eigenvalues precisely coincide with the roots of a 
certain system of algebraic equations, which arises as a part of the coordinate [2] or algebraic [3] 
Bethe Ansatz. 

The underlying algebraic structure behind the construction of the commuting transfer ma- 
trices, termed T-matrices or T-operators, is by now well understood. It is connected with the 
simplest representations of quantum groups, which are closely related to the standard finite- 
dimensional representations of classical Lie algebras. In contrast, the algebraic construction of 
the Q-operators appears to be a more complicated and, at the same time, more interesting prob- 
lem. Much progress in this direction has already been achieved from a case-by-case study of 
various models, see e.g. [H[2l|4l[5j|6l[7l|8j|9l [TUl IH] but the problem still continues to reveal its 
new features. 

In this paper we develop a new approach to Q-operators by connecting them to the theory of 
Yangians, which are the simplest examples of quantum groups. In doing so we shall develop new 
aspects of the theory of infinite-dimensional representations of Yangians, naturally leading to a 
systematic and transparent construction of the Q-operators. Here we illustrate our approach on 
the compact g[(n)-spin chains, but the results may be readily generalized to other models, and 
in particular to supersymmetric spin chains [12J. 

Let us then consider the integrable g[(n)-spin chain with the well-known Hamiltonian 



in the presence of twisted periodic boundary conditions, 



where <E>i, <3? 2 , • • • , $ 

n is a set of fixed twist parameters (or fields). Here e a j, denotes the n x n 
matrix unit (e a &)jj = 5 a i5bj and the superscript "(£)" refers to the quantum space of the l-th. 
spin in the chain. Each "spin" can take n different values a = 1,2, ... , n. It is easy to check 
that for the quasi-periodic boundary conditions (| 1 . 2 [) the numbers m\,mi, . . . , m n , counting the 
total number of spins of type "1", "2", . . . , "n" in the chain, are conserved quantum numbers 
for the Hamiltonian (II. ip . Due to these conservation properties this integrable model can be 
solved via the "nested" Bethe Ansatz technique |13[ [14] , which leads to the well known result for 
the eigenvalues of (jl.ip . see (I6.8f> below. They are expressed through solutions of the already 
mentioned algebraic equations, commonly called Bethe Ansatz equations. 

It is relatively well known that there are different but equivalent forms of the Bethe Ansatz. 
In fact, it is easy to argue that there are precisely n\ different Bethe Ansdtze in our case, related 
by all possible permutations of the occupation numbers 7714,777-2, • • • ,777 n . Indeed, there are n 
ways to choose the bare vacuum state, then n — 1 ways to proceed on the second "nested" stage 
of the Bethe Ansatz and so orQ. These options can be conveniently depicted by directed paths 

x One should keep in mind, however, that the above argument only applies to the case where all fields $ a take 





2 



on a Hasse diagram which spans an n-dimensional hypercube. The nodes of the hypercube are 
labeled by increasing integer sets / = {ai, 02, . . . , a p } C {1, 2, . . . , n}, where < p < n. See the 
n = 3 example in Fig. [TJ-a. There are exactly 2 n nodes on the diagram and exactly n! ordered 
paths from the bottom to the top. Then each path is related to a particular variant of the 
Bethe Ansatz, while nodes on that path are related to the so-called Q-functions entering the 
corresponding Bethe Ansatz equations. This concise description was proposed in |16j . Note that 
the usual nested Bethe ansatz proceeds from the top to the bottom of this diagram. 



Q{1,2,3} 




Figure 1: a) Hasse diagram for 0t(3), b) An elementary quadrilateral, whose bottom node corre- 
sponds to the set / = {ai , 0,2, ■ ■ ■ , a p }. 

Our approach precisely reproduces this picture. We explicitly construct the 2 n different 
Q-operators, corresponding to the nodes of the Hasse diagram, and subsequently derive the 
nested Bethe Ansatz equations and eigenvalues of the Hamiltonian without any reference 

to the eigenvector construction. The different Q-operators are not functionally independent. 
They satisfy the so-called Hirota equations, defined on the direct product of the n-dimensional 
hypercube and the real line. The equations have the same form for every quadrilateral of the 
Hasse diagram (see Fig. [IJb), 

&{a,b}QlUaUb(z) Ql(z) = QlUa{z ~ \) QlUb{z + \) ~ QlUb(z ~ \) Q/Ua(> + §), (1-3) 

where 

A {ajfe} = 2i sin ( (1.4) 

and the expression I U a denote the union of the sets / and the one-element set {a}. The 
reader might be aware that the Hirota equations frequently arise in the analysis of quantum 
integrable models, see e.g. |17] , The equations (11 .3j> for the Q-operators and their eigenvalues 
were gradually developed for various models related to gl(n) algebra in (6J [181 E3 120] . The form 

generic, non-zero values. If some or all of the fields vanish, or else take certain special values, only a few of the 
Bethe Ansatze are well defined, while other ones typically suffer from multiple roots or the so-called "beyond the 
equator" problem |15j . 
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of the operatorial equations (jl,3p presented here exactly coincide with that of the eigenvalue 
equations of ref. [16p. 

The most essential algebraic property of the Q-operators is concisely expressed by the single 
factorization relation 

T+(z, A n ) ~ Q {1} (z + X[) Q {2} (z + A' 2 ) • • • Q {n} (z + X'J (1.5) 

where T + (z,A n ) is the transfer matrix associated with a highest weight infinite-dimensional 
representation (Verma module) with the $j[(n)-weights A n = (Ai, A2, • • • , A n ). Here = A n + 
p n denotes the weights shifted by the half sum of the positive roots, p n of the algebra sl(n). 
Remarkably, this relation alone allows to derive all functional relations, satisfied by various 
"fusion" transfer matrices and Q-operators. As a result the model may then be solved in an 
entirely algebraic fashion. In particular, the transfer matrices 

T(z,A n ) =*det||Qi(z + A$)||i<ij<n (1.6) 

associated with finite dimensional representations, are expressed as the determinant of a matrix 
constructed from Q-operators. This approach was originally developed mainly for field theory 
models connected with quantized (or g-deformed) affine algebras [6j [22], [19l [23| 124"] . 

In |llj we applied the same idea to solve the famous Heisenberg XXX spin chain corresponding 
to n = 2 in (11. ip . which is much simpler than its g-deformed counterparts. This allowed us to 
observe some new algebraic structures which were obscured or did not manifest themselves for 
field theory models. We found that for the XXX model the relation (ll.5|) arises from very elegant 
factorization properties of the st(2)-invariant Lax operator 

£+(*, A 2 ) = z + ^, a »J»=( Z +/ J " j3 ) , (1.7) 

a=l ^ ' 

which acts on the tensor product of a spin-^ module C 2 and the infinite-dimensional highest 
weight representation of sl(2) with an arbitrary spin j. The Sl(2) generators J 3 , 3^ in this case 
are realized in the Holstein-Primakoff form, 

J- =b f (b t b-2j) , J+ = -b, J^j-b+b, (1.8) 

where b^ and b are generators of the oscillator algebra 

H: [b ,b f ] = 1, h = b f b +i . (1.9) 

Here we use Ql(2) representation labels, such that A2 = {j, —j)- The factorization in question 
involves the two simpler "constituent" L-operators 

with two sets of mutually commuting oscillators. The precise factorization formula reads 

2 There are also various sypersymmetric extensions of (|1.3p . the related bibliography can be found in [161 121] , 
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where 



z\ 



z + j + 



Z2 



1 

2' 



and the generators J 3 , J 1 * 1 , are realized in the Holstein-Primakoff form 
replaced by b{ and bi. It is useful to rewrite (jl.lip in a compact form 



Li(zi)L 2 (z 2 ) =S£ + (z,A 2 ) 



7-1 



where 



b 2 
1 



(1.12) 
with b^ and b 

(1.13) 
(1.14) 



A similar formula for the reversed order product L 2 Li can also be found in |llj . These factor- 
ization equations were used to derive the n = 2 version of (|1.5p as well as the complete hierarchy 
of all functional relations, leading to a new algebraic solution of the XXX magnetic spin chain. 

It is interesting to note that the LHS of (jl.lip contains two operators (jl.lOp . which are first 
order polynomials in the spectral parameter z, and, at the same time, their product in the RHS 
is also a first order polynomial in z. This is explained by the fact that in (jl.lOp terms linear in 
z are proportional to two different degenerate matrices whose product is equal to zero. Thus, in 
contrast to (|1.7p . where the term containing z is proportional to the unit matrix, the operators 
(jl.lOp start from degenerate matrices. 

By extending this key observation to the g((n) case, we have completely classified all first order 
L-operators and studied their fusion and factorization properties. In general they correspond 
to some special infinite-dimensional representations of the Yangian Y(qI(ti)) related to a direct 
product of gl(p) ® 7^P( n_ P) ) involving the algebra gl(p) with p < n, and a number oscillator 
algebras %. In particular, we find that the sl(n) Lax operator 



£ + (z,An) 



z + ^egb 

a.b 



>ba 



(1.15) 



for infinite-dimensional highest weight representations of the flt(n) generators J a b factorizes into 
n partonic Lax operators L a (z) (a = 1, . . . , n), 



( 1 



b 



a.l 



a, a— 1 



»l,a 



— "a-l,a 
z -h a b 



\ 



b 



a+l,o 



t 

a,a+l 
1 



V 



b f 



1 / 



;i.i6) 



which is the generalization of (jl.lOp for the gl(n) case. Each of these operators contains n — 1 
independent oscillator pairs (b& , bL), b = 1, 2, . . . , n, b ^ a, 



b^a 



The factorization formula, generalizing (jl,13p . now reads 

U{z + A;)L 2 (z + A' 2 ) • • • L n (z + A^) = S c C + (z I An) G C S C \ 



;i.i7) 



:i.i8) 
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Here the shifted weights A' n = A n + p n are as in (jl.5p , p n is given by f|4. 18j) and the quantities Sc 
and Gc are generalizations of those in (|1.14p . Eq. (jl,18p immediately implies the factorization 
relation (jl,5p . 

Finally, note that the formula (|1.16p (and its generalization (|2.20p ) define extremely simple 
first order solutions of the Yang-Baxter equation. The solutions are new and were not previously 
considered in the literature. We will now proceed to derive and precise all above statements. 

2 The Yang-Baxter Equation and Representations of Yangians 

The Hamiltonian (II. ip with twisted boundary conditions commutes with a large commuting 
family of T- and Q-operators. In this paper we explicitly construct these operators via traces of 
certain monodromy matrices associated with infinite-dimensional representations of the harmonic 
oscillator algebra. To do this we need to find appropriate solutions of the Yang-Baxter equation 

R(*i - z 2 ) (L(zi) ® l) (l <8> L(z 2 )) = (l <g> L(z 2 )) (h(zi) ® l) R(*i - z 2 ), (2.1) 

where R(z) is an n 2 x n 2 matrix, 

R(z) : C n ® C" -> C n <g> C n , R(z) = z + P , (2.2) 

acting in the direct product of two n-dimensional spaces C n ® C n . The operator P permutes 
the factors in this product. The operator L(z) is an n x n matrix, acting in a single copy of 
the space C n , whose matrix elements are operator- valued functions of the variable z belonging 
to some associative algebra y. To be more precise, the Yang-Baxter equation (|2.1D provides 
defining relations of the Yangian algebra y = Y(Ql(n)), introduced by Drinfeld [25] (for a recent 
comprehensive review see |26j). 

Let Lij(z), i,j = 1,2, ... ,n denote the matrix elements of L(z). From ([2.ip it follows that 

(y - x) [Lij(x), L ki (y)] = L kj (x) Lu(y) - L kj (y) Ltf(x), i, j, k, £ = 1, 2, . . . , n. (2.3) 

Writing as a Laurent series with operator- valued coefficients 

L lJ (z)=L^ + Lffz^ + Lffz- 2 + ..., (2.4) 
one obtains an infinite set of commutation relations 

min(r,s) 

[4;ua= e {Ltr a)L u~ l) -< rX) ^ +s ~ a) ) ( 2 - 5 ) 

0=1 

( r \ 

for the elements L\- , r = 0, 1, 2, ... , oo. Thus the problem of solving Eq. (|2.ip reduces to the 

construction of representations of the infinite-dimensional quadratic algebra (12. 5p . In our ap- 
proach the most important role is played by the simplest representations, where the series (12. 4p 
truncates after the second term. To within a trivial change in the normalization these repre- 
sentations correspond to L-operators, which are first order polynomials in the parameter z. All 
previously known L-operators of this type can be brought to the form 

Cij(z) = z 5ij + Jji , (2.6) 

where Jij, i, j = 1, 2, . . . , n denotes the standard set of generating elements of the algebra gl(n), 

[Jij,Jki] = bkjJil—bitJkj- (2.7) 



G 



The Yang-Baxter equation for the L-operator (|2.6p is satisfied on the algebraic level by virtue of 
the commutation relations (|2.7p . Therefore one can choose in (|2.6j) arbitrary gl(n) representations 
for the generators Jij. It is obvious that the addition of a constant to the spectral parameter z 
in (|2.6p can be compensated by the subtraction of the same constant from the central element 

C n = JlX + J22 + • • • + Jnn (2-8) 

of the algebra (I2.TI) . Therefore it is tempting to eliminate this spurious degree of freedom by, for 
instance, imposing the condition C n = and restricting (|2.7p to the algebra sl(n). Here we will 
not do so, but will instead work with the full algebra fll(n). This is helpful for clearly exposing 
the Weyl group symmetry of the problem at hand. 

We will now show that, excitingly, there exist further first order L-operators, different from 
(|2.6p . We will present their complete classification. To begin, let us recall a simple symmetry of 
the Yang-Baxter equation (12. ip . The i?-matrix (|2,2p is GL(n)-invariant in the sense 

R(z) = (G®G)R(z) (G&G)- 1 , G G GL(n), (2.9) 

where G is any non-degenerate nxn matrix. It follows that if L(z) satisfies (I2.ip . then any other 
operator of the form 

L(z)=FL(z)G, F, G G GL(n), (2.10) 

with arbitrary F,G£ GL(n) satisfies again the same equation. Furthermore, the matrices F, G 
may contain operator- valued matrix elements, as long as these commute among themselves and 
with all other elements of L(z). 

From (12. 5h it then immediately follows that the elements 1$ are central, i.e. they commute 



among themselves and with all for r > 1. Therefore, we may reg ard L(°) numerical 
nxn matrix. Applying the transformations f|2. 10|) . this matrix can always be brought to diagonal 
form 

L® =diag( 1,1, ...,1,0,0, ...,0), p = l,2,...,n, (2.11) 

p-times (n— p)-times 

where p is an integer 1 < p < n. The number p coincides with the rank of the matrix . It is 
invariant under the linear transformations (|2,10p ). Evidently, if p = n, the leading term in the 
series expansion (|2.4p is the unit matrix. This case is well studied in the existing representation 
theory. In fact, the assumption that the series (12 . 4|) starts with the unit matrix is usually 
included into the definition of the Yangian. Here we will not make this assumption, and will 
consider instead the more general case with arbitrary 1 < p < n. 

Let us concentrate on the simple case when the series (|2,4p truncates after the second term, 



(r) 

i.e. assume that all L\- = for r > 2. It is convenient to write the only remaining non-trivial 
coefficient L^ l > block matrix 



Kb \ B ah 



Cab \ D ah 



(2.12) 



where A, B, C and D are operator- valued matrices of dimensions p x p, p x (n — p), (n — p) x p 
and (n — p) x (n — p), respectively. We furthermore assume that all undotted indices run over 
the values {1, 2, . . . ,p}, whereas their dotted counterparts take on the values {p + 1, . . . , n}. 

l<a,b<p, p+l<a,b<n. (2.13) 
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Substituting (|2.1ip and (I2.12p into (|2.5p . one immediately realizes that the elements D.i are 
central, i.e. they commute among themselves and with all other elements of ZA 1 ). The other 
commutation relations read 



[ B c&> C db] 



&ae A c b 5 c b A ae , 
Sab D hh , 



\A a bi Bcc\ 
l B ab' B ce] 



-he B a 



0, 



[A a b , C c , 

[C&6; C c , 



. 



(2.14) 



Using the remaining freedom of making transformations (|2.10p , which do not affect the form of 
IA°) in ([2.1ip . one can then bring the matrix D to diagonal form with zeroes and ones on the 
diagonal, in similarity to (|2.1ip . Here we are only interested in highest weight representations 
of the algebra (|2.14p . These representations admit a definition of the trace, as required for the 
construction of transfer matrices in Section H] below. For this reason we only need to consider the 
non-degenerate cased, det where the diagonal form of D coincides with the (n— p) x (n—p) 

unit matrix 

D ab = S ab> P+l<d,b<n. (2.15) 

The resulting algebra (|2.14p can be realized as a direct product of the algebra Ql(p) with p{n—p) 
copies of the harmonic oscillator algebra: 



(2.16) 

p + 1, ... ,n and b = 
(2.17) 



Introduce p(n — p) independent oscillator pairs (bit, bL), where i 
I, . . . ,p, satisfying the relations 

KbMel = SdeSbc- 

Furthermore, let J a b, a,b = 1,2, ... ,p denote the generators of the algebra gl(p) defined by (12 .7p . 
where n is replaced by p. The generators J a \, commute with all oscillator operators in (|2,17p . The 
connection of the algebra (12. 14|) with the product (I2.16P is established by the following relations 



A 



ab 



J, 



ab 



(bLb a + W, B ab = *>l b i Cab = -b ab , (2.18) 
where the upper bar in the notation J ab denotes the transposition of the indices a and b, 



E 

c=p+l 



iJ ab — <Jba- 

The corresponding L-operator can be written as a block matrix, 

/ 



(2.19) 



L {l,2,...,p}(^) 



+ J ab 



n 

E 

c=p+l 



(bLb cfe + 



V 



ab 



b 



ab 



ab 



(2.20) 



/ 



where the rows are labeled by the indices a or a and the columns by b or b, in similarity to 
(|2,12p . Note that the p x p matrix of the generators of which enters the upper left block, 

is transposed, i.e. the a-th row and 6-th column in this block contains the element J ab = J\, a - 



3 It appears that for det D = the algebra (|2.14[) does not admit a definition for a suitable trace as needed for 
the construction of transfer matrices commuting with the Hamiltonian (|l.ip . 
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The matrix (|2.20p contains the parameter z only in its first p diagonal elements. By simul- 
taneous permutations of rows and columns in (|2.20p one can move these z-containing elements 
to p arbitrary positions on the diagonal, labeled by a set of integers I = {a±,a2, ■ ■ ■ ,a p }. We 
shall denote the L-operator obtained in this way by Jj j(z). Within this convention the operator 
(|2,20p corresponds to the set I = {1, 2, . . . ,p}, as indicated by the subscript in the LHS of this 
equation. 

The "partonic" L-operator (j 1 . 161) is a particular case of (I2.20P with p = 1, while the standard 
L-operator (12. 6h corresponds to p = n. We would like to stress that for p < n the formula (I2.20P 
yields novel solutions of the Yang-Baxter equation (|2.ip . The only exception is the simple case 
n = 2, p = 1, where this solution was previouly known \27\ [28| [29]). Note, also that the n = 3 
solutions can be obtained in the rational limit of trigonometric solutions obtained in [301 131] . 

The formula (I2.20P provides an evaluation homomorphism of the infinite-dimensional Yangian 
algebra (I2.5P into the finite-dimensional algebra (|2.16p . 

Y(gl(n)) -> fll(p) ® H®p("-p) , l<p<n. (2.21) 

This means that for any representation of this finite-dimensional algebra the equation (I2.20P 
automatically defines a representation of the Yangian and a matrix solution of the Yang-Baxter 
equation (|2.ip . Conversely, any first order matrix L-operator with a rank p leading term 
and a non-degenerate matrix D in (I2.12P is, up to a transformation (I2.10p . equivalent to the 
canonical L-operator (|2.20p with some particular representation of the algebra (|2.16p . It is worth 
noting that the transformation (I2.10P 

L { i,2,..., P }(*) -»• FL {1)2i ... iP} (;z)F- 1 , (2.22) 

where F is a, block diagonal ma,trix containing the matrices Wp G GL^p^j, W n —p G GL(n - p) on 
the diagonal, leaves the form of (|2,20p unchanged. 

In preparation for a necessary analysis below we need to introduce some notation for the 
irreducible highest weight representations of gl(n). The highest weight vector vq is defined by 
the conditions 

Ja,a+iv = 0, a = l,2,...,n-l. (2.23) 
Here we will use g[(n)-type representation labels^ A n = (Ai, A2, • • • , A n ) where 

JaaV = \ a vo, a = l,2,...,n. (2.24) 

We will denote by ir^ the corresponding infinite-dimensional highest weight representation with 
arbitrary weights, and by ir\ n the finite-dimensional irreducible representation with 

Ai>A 2 >...A n , A a — At G Z . (2.25) 

The analysis of this section extends the previous results of |27] devoted to n = 2 case. The 
properties of the finite dimensional representation of the Yangian Y(g[(2)) associated with the 
L-operator (12. 6p can be found in |32[ 133"! 126] . 

4 The s((n)-type weights fi = (jui, . . . , ^n-i), 

fj. a — A a — A0+1, a = 1, 2, . . . , n — 1 , 

are inconvenient from the point of view of the Weyl group symmetry. They will not be used in this paper. 
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3 Fusion and Factorization of L-operators 



An essential part of our analysis in the following is based on some remarkable decomposition 
properties of the product of two L-operators of the form (|2.20p . The Yangian y = Y(gl(n)) is a 
Hopf algebra, see e.g. |26| . In particular, its co- multiplication 

y^y®y (3.1) 

is generated by the matrix product of two L-operators, corresponding to two different copies of 
y appearing on the RHS of (I3.1|) . 

Our main observation is related to the co- product of two operators L/(z) and ~Lj(z), defined 
by (I2.20p for two non-intersecting sets I D J = 0, 

L(z) = L [ j ] (z + zi)L [ j ] (z + z 2 ), (3.2) 

where the quantities z\ i denote arbitrary constants. Let the sets I and J contain p\ and p 2 
elements, respectively. It is obvious that 

(i) the product (|3.2p is of the first order in the variable z, and that 

(ii) the matrix rank of the term linear in z in fj3.2j) is equal to p\ + P2- 

The simple meaning of the co-multiplication is that the matrix product of two L-operators, each 
of which satisfies by itself the Yang-Baxter equation (|2.ip . solves this equation as well. All 
solutions which possess the above properties (i) and (ii) were classified in the previous section. 
Therefore, by using a transformation of type (|2.9p , the RHS of (|3.2p can be brought to a particular 
case of the canonical form (I2.20P with p = p\ + pi- It turns out, however, that the expressions 
for the matrix elements of the resulting L-operator are rather complicated and their explicit 
connection to those of (I2.20|) is far from obvious, even though these elements satisfy the same 
commutation relations. In order to make this connection more transparent we apply a suitable 
operatorial similarity transformation S to each matrix element such that it rearranges the basis 
of the oscillator algebras contained in An pi <S> An P2 ■ Furthermore, the formula (I3.2p contains 
two constants z\ and Zi- Obviously only their difference is an essential parameter, whereas the 
sum may be absorbed into the spectral parameter z. Therefore, without loss of generality, one 
can set 

zi = A + f, z 2 = -f, (3.3) 

where A is arbitrary. This particular parametrization is chosen to simplify the subsequent for- 
mulae. 

Proceeding as described above, one obtains, 

L(z) = L™ (z + A + f ) \}f (z-^)=S (L 7uJ (z) G)s-\ (3.4) 

where G is a z-independent matrix, whose elements commute among themselves and with all 
elements of ~L Iu j(z). It should be stressed that the resulting L-operator ~L I{J j(z) is only a special 
case of (I2.20p . since it is connected to some specific realization of the algebra A n ,p 1 +p 2 m terms 
of the direct product of the two algebras A ntPl ® A n ^ P2 as defined in (|2.16p . Note that the RHS 
of (I3.4p is of course a particular case of the transformation (|2,10p with F = 1. The explicit 
expression for the matrices appearing in (|3 . 4[) are presented below. Some additional details of 
calculations are given in the Appendix lAl 
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By permuting rows and columns any two non-intersecting sets I and J can be reduced to the 
case when I = {1, . . . ,pi} and J = {p\ + 1, . . . , p\ + p2\. So it is suffficient to consider this case 
only. Introduce three types of indices 



a,b,el, a,b e J, a, b G {pi +P2 + 1, • 
It will be convenient to rewrite (|2,2U|) as a 3 x 3 block matrix 

/ 



,n} 



(3.5) 



J {l,2,...,pi} 



^a 6 + j2-E( b - ] b 2 +5**) 


b t[l] 
ao 


b t[l] 

ab 








ao 







"db 





^do 



, (3.6) 



/ 



where the size of the diagonal blocks is equal to p\ x pi, p2 x j>2 and (n — p\ — p2) x (n—pi — P2), 
respectively. The superscript " [1]" indicates that the corresponding operators belong to the "first" 
algebra in the co-multiplication (|3.ip , which in the considered case is realized, by the algebra 1An.p1 



[21 

write Lj (z) as 






$ab \ 







D db ; 2 °db + ^ dfe 




b tP] 

do 


! 


do 


5 db 



I 



(3-7) 

where superscript "[2]" labels operators from the "second" algebra, which is the algebra A ntP2 . 
By construction, all operators labeled by the superscript "[1]" commute with those labeled by 
the superscript "[2]". Recall also that the algebra (|2.16p has a direct product structure, so the 
generators j\ \r} and J^ 2 } commute with all oscillator operators. 

a0 db 



With the notation introduced above the similarity transformation S has the form 

S = S1S2 , 



where 



and 



5i=exp(EE b ^ b ii 21 )' 



ft— plEEE^bJM 

cGJ cS J C0U.J 



(3.8) 
(3.9) 

(3.10) 
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The matrix G has the form 



V 



8ab 









ab 





6 db 














(3.11) 



/ 



The similarity transform S\ serves to expose the fact that the matrix entries of L/uj(z) commute 
with the entries of G. The similarity transform £2 brings hnjj(z) to the form (|2.20p . 

Finally, we want to write the operator ~L I{J j(z) in (|3,4p in the form (|2.20p with p = pi +P2- To 
do this we need to make the following identifications for the generators Jij, i,j = l,...,pi+p2 
in the upper diagonal block of (|2.20p 



t [1] uW 



CgJ 



J db- J db + Z^°cd ° bc ' 



t [1] uW 



eel 



(3.12) 



ab 



tab 



\cEl c&J / 



Ab 



t[i] 

ab 



+E^3"S. w -E b Ii ,| 7 



t [i] -Tf [i] 



eg J 



eg/ 



where we have used the convention (|3.5p for numerating indices and Jjj = Jjj. 

Furthermore, let the indices A, B run over the values 1, 2, . . . ,pi + P2 and A, B over the values 
Pi + P2 + 1) • • • j n - Introduce operators 



b [1] 


Bel, 


AB ' 


b [2] 


B G J, 


AB ' 


from 


([3~H) can 



'ab ' 



A £ I, 



b[f, A £ J. 



(3.13) 



/ Z5 AB + J; 



J {l,2,...,pi+p 2 }( Z ) 



E 

C0IUJ 



- C ac C cb + 2^ AB ) 



V 



J 



(3.14) 



which has the required form as in (|2.20p . 

The formulae f|3. 1 2[) give a homomorphism of the algebra Ql(pi + P2) into the direct product 



0l(pi + P2) ->- 0l(pi) ® fll(p2) ® U® P1P2 = B. 



P1,P2 1 



(3.15) 



which for p\^, 7^ has only infinite-dimensional representations (similar representations appeared 
in |34j). An important feature of this map is that if one choses highest weight representations for 
both algebras 0t(pi) and Qi{p2) then the formulae (|3.12p define a highest weight representation 
of Ql{pi +P2)- It is easy to check that the conditions (|2.23p are satisfied on the product of the 
corresponding highest weight vectors Vq \ and the standard Fock vacuum for all oscillator 
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algebras appearing in (|3.12p . The Ql(p± + p 2 )-weight of the resulting representation is easy to 
obtain from f|3. 12[) 

A P1+P2 = (a[ 1] + A, A 2 1] + A, . . . , AW + A, \f\ A 2 2] , . • • , Ag) , (3.16) 

where A is an arbitrary parameter from (|3.4I) . 

The L-operators in the first product in (I3.4|) have the superscripts [1] and [2], which indicate 
that they belong to different algebras (|2,16p with p = pi and p = P2, respectively. By the same 
reason it is useful to rewrite the RHS of (I3.4j) supplying similar superscripts 

<m=s{L [ SA z ) Gl2 ' ] ) s ~\ ( 3 - 17 ) 

where the superscript [1 ; ] indicates the algebra f|3. 15|) and the superscript [2'] indicates the 
product of oscillator algebras %®P!P 2 . Note that the matrix G could be considered as a z- 
independent L-operator, also satisfying the Yang-Baxter equation (|2.ip . In view of this Eq. (|3.17p 
also describes the co-multiplication of two representations of the Yangian. 

Consider now some particular consequences of formula (|3.4p . Using it iteratively with p\ = \ 
and taking into account (I3.16P one can obtain an arbitrary product of the elementary L-operators 
(|1,16|) . Let I = (ox, . . . , a p ) be an ordered integer set, 1 < a\ < . . . < a p < n, and Cf(z | A p ) a 
specialization of the L-operator ()2.20[) to the infinite dimensional highest weight representation 
7r^ of the algebra gl(p) 

L+(z|A p ) =tt+ [Lj(z)], Ap = (X 1 ,X 2 ,...,X P )- (3-18) 
Define also the shifted weights (cf. (|4.18p ) 

a; = (a , 1 ,a' 2 ,...,a;), x'. = \ j + ^M±i y j = i,..., P . (3.19) 

Then it follows that from (I3.4p 

L ai (z + A;)L a2 (z + A' 2 ) • • • L ap (z + A;) = Sj Lj(z \ A p ) Gj5f 1 (3.20) 

where the matrices Sj and Gj are products of the expressions of the type (|3.8p and ([3.1ip arising 
from the repeated use of the formula (|3.4p . In the particular case p = n the last formula provides 
the factorization for the L-operator (|2.6D . 

C + (z\A n ) = 7r+ n [£(z)], (3.21) 

evaluated for the infinite-dimensional highest weight representation tt^ in the auxiliary space, 

Li(* + A;)L 2 (z + A' 2 ) • • • L n (z + A;) = S c C + {z \ A n ) G C S C \ (3.22) 

An independent proof of this fact is given in the appendix |Bj 

4 Construction of the Q-operators 

The purpose of this section is to define the T- and Q-operators. They have to commute with 
the Hamiltonian (jl.ip of the twisted compact g[(n)-spin chain of length L. These operators act 
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on the quantum space V, which is an L-fold tensor product of the fundamental representations 
of the algebra Ql(n), 

V = C n ®C n ®---®C n . (4.1) 



L— times 



As before, solutions of the Yang-Baxter equation (I2.1h are considered as n by n matrices, acting 
in the quantum space of a single spin. Their matrix elements are operators in some representation 
space V of the Yangian algebra (I2.3|) . This representation space will be called here the auxiliary 
space. For each solution of (|2.1|) one can define a transfer matrix, 

T v (z) =Tiv|BL(*)®L(3)®...®L(z)}, (4.2) 

where the tensor product is taken with respect to the quantum spaces C n , while the operator 
product and the trace is taken with respect to the auxiliary space V. The quantity D is a 
"boundary twist" operator acting only in auxiliary space, i.e. it acts trivially in the quantum 
space. This boundary operator is completely determined by the requirement of commutativity 
of the transfer matrix f|4. 2|) with the Hamiltonian (|1.1|) , which leads to the following conditions 

D (Mz)) ab B- 1 = e l (*>-*»> (Hz)) ab , a,b = l,...,n, (4.3) 

where the fields $ a are the same as in (|1.2p . Note that these fields enter (|4.3p only through their 
differences. It is convenient to set 

n 

X>a = 0. (4.4) 

0=1 

Solving (I4.3|) for the general L-operator (|2,20j) with an arbitrary set I = {a%, a^-, ■ ■ ■ , a p }, one 
obtains 

Vj = exp {* £ d> a J aa d> fc ) bl b b ba \ . (4.5) 

In the following we will use some important properties of the trace over the Fock representa- 
tions of the oscillator algebra 

[b,b t ] = l, [h,b] = -b, [11,^]=^, h = btb + i. (4.6) 

This algebra has two Fock representations, 

T+ : b|0)=0, |fc + 1} = b f (4.7) 

and 

JF_ : b t |0) = 0, \k + l) =b\k), (4.8) 

spanned on the vectors \k), k = 0, 1, 2, . . . oo. These representations can be obtained from each 
other via a simple automorphism of (14.6|) . 



b^-b f , b f ^b, h^-h. (4.9) 

Let P(b,bt) be an arbitrary polynomial of the operators b and b^. Below it will be convenient 
to use a normalized trace over the representations J r ±, 

Tr^{e i$h P(b,bt)} 

TrJe^ h P(b,b^)} d = -, J r = J r ±, (4.10) 
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where T is either or T-, and Trjr denotes the standard trace. An important feature of 
the normalized trace (|4,10p is that it is completely determined by the commutation relations 
(14, 6p and the cyclic property of the trace. It is therefore independent of a particular choice 
of representation as long as the traces in the RHS of (|4.10p converge. Alternatively, one can 
reproduce the same result by using explicit expressions for the matrix elements of the oscillator 
operators in (|4,10p . Then the trace over J- + converges when Im<3> > and the trace over T— 
when Im $ < 0. Both ways of calculation lead to the same analytic expression for the normalized 
trace. Thus it is not necessary to specify which of the two representations T± is used. 

We are now ready to define various transfer matrices all commuting with the Hamiltonian 
(jl.ip . Consider the most general L-operator (|2.20p with an arbitrary set / = {ai,ci2, ■ ■ ■ ,a p }, 
where p = 1,2, ... ,n. Recall that the matrix elements of (|2,20p belong to the direct product 
(I2.16P of the algebra Ql(p) and of p(n —p) oscillator algebras. Choose a finite-dimensional repre- 
sentation 7T A with the highest weight A p for the fj[(p)-factor of this product. Then substituting 
(I2.20P and (14.5P into (14. 2p one can define rather general transfer matrices 

X T (z, Ap) = e te C£. e /*«) l^T*^ (n _ ri {Mj(*)} , (4.11) 

where M/(z) is the corresponding monodromy matrix, 

M T (z) = Vj L/(z) <g> Li(z) <S> ■ ■ • <g> Li(z) . (4.12) 

Here Tr j-p(n-p) denotes the normalized trace (I4.10P for all involved oscillator algebrasEI, while 
Tr^ denotes the standard trace over the representation ir\ p of $l(p)- The exponential scalar 
factor in front of the trace is introduced for later convenience. 

Similarly one can define a related quantity where the g[(p)-trace is taken over an infinite- 
dimensional highest weight representation 7T^ , 

X+(z, A„) = e "(£.€x*»> { Mj (^)} , (4.13) 

while the rest of the expression remains the same as in (14.1ip . Note that in the case of (14.1ip the 
weights A p = (Ax, A2, . . . , X p ) satisfy the conditions (I2.25p . In contradistinction, in (|4.13p these 
weights are arbitrary. 

In the limiting case p = n, the general L-operator (|2,20p simplifies to (|2.6p . while the expres- 
sion ([4.5p simplifies to 

n 

V = ^{1,2,...,„} = exp {» ^aJaa) ■ (4.14) 

o=l 

In this case the definition (|4.1ip reduces to that for the standard T-operator 

T Ab (z) = X {h2 _ n} (z,A n ) = Ti VAn \vC{z)®C{z)---®C{z)}, (4.15) 

associated with the finite-dimensional representation 7TA n of the algebra 0[(n) in the auxiliary 
space. Here C(z) denotes the L-operator (|2.6p . Likewise, the formula (I4.13P reduces to the 
T-operator 

T A>)= X {l,2,...,n}(^ A n) ( 4 - 16 ) 

5 It is easy to check that all possible expressions under the trace in (|4.1ip for each oscillator algebra are exactly 
as in the LHS of (|4.10[) for some polynomial P and some value of 4>. Thus the definition (|4.10|l is sufficient to 
calculate all oscillator traces in (|4.11[) . 
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associated with the infinite-dimensional representation 7r^ . The above two T-operators are con- 
nected due the Bernstein Gel'fand Gel'fand (BGG) resolution of the finite dimensional modules 
|35| . The BGG result allows one to express finite dimensional highest weight modules in terms 
of an alternating sum of infinite dimensional highest weight modules. This implies that the 
T-operator (|4.15p for a finite dimensional module can be written in terms of ([4.16P as 

a&Sn 

where p n is a constant n-component vector 

Pn= (V^,---, 1 ^), ( 4 - 18 ) 

which coincides with the half sum of the positive roots of the algebra sl(n). The summation 
in (|4.17p is taken over all permutations of n elements, a £ S n , and 1(a) is the parity of the 
permutation a. The relation (|4,17p and its connection to the BGG resolution were first obtained 
in |19j in the context of U q (sl(3)), while the n = 2 case was previously considered in [6j I22j . 
Similarly, for ([4. lip one has 

X I (z,A p )= £(-l)^>X+(z ) a(A p + p p )- j0p ), (4.19) 

where p p is a p-component vector defined as in f|4. 18[) with n replaced by p. 

Another limiting case of (I4.1ip corresponds to the representation tt^ turning into the trivial 
one-dimensional representation of Ql(p) with weight A = (0, 0, ... , 0). As we shall see below the 
resulting operators 

Qj(z) = X T (z, (0)) (4.20) 

are actually the Q-operators, whose eigenvalues appear in the nested Bethe Ansatz equations. 
Let us enumerate these Q-operators. It is convenient to start formally from the exceptional case 
p = 0, corresponding to an empty set 1=0. By definition we set 

Q (z) = 1. (4.21) 

For the next level p = 1 there are obviously n sets / consisting of just one element / = {a}, 
a=l,2,...,n. The general L-operator (|2,20p in this case takes the simple form f 1 1 . 1 6 [) and the 
twist operator (|4.5|) simplifies to 

V a = V {a} =exv{-iY^($a-$c)bl 6 b 6a }, a = l,2...,n. (4.22) 

In this way one obtains from (|4.1ip 

Qaiz) = Q {a} (z)=X {a} (z,{0))=e iz ^ Tf^„_i ) {D L (2;)®L (2)(8)---®L a («)} > (4.23) 

where a = 1,2 ... ,n and L a (z) is given by ([1.16p . 

More generally, for the level p there are (") increasing integer sets / = {ai,...,a p } C 
{1, 2, . . . , n}, which numerate the Q-operators (|4,20p . For the highest level p = n the definitions 
(I4.20P and (14.1ip immediately lead to the result 

Q{i,2,...M z ) = zL > ( 4 - 24 ) 
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where L is the length of the chain. Altogether there are 2™ different Q-operatora^| , including 
(|4.2ip and (|4.24p . These operators form in conjunction with all T and X-operators a commuting 
family and therefore can be simultaneously diagonalized. It is easy to see that their eigenvalues 
have the form 

rrij 

Q J (z)=e"(E o6 ,*.) H(z-zl), mi = ^2m a , (4.25) 

k=l a£l 
where, for each eigenstate, the numbers m a are the conserved occupation numbers, 

mi + m 2 + • • • + m n = L , (4.26) 

defined in the introduction. 

We would like to stress, that in general the operators ~X.j(z,A p ), defined in (|4,lip . involve 
the trace over a representation of the Lie algebra Ql(p) and the trace over a number of Fock 
representations of the oscillator algebra, whereas the T-operators involve only the Lie algebra 
trace and the Q-operator only the oscillator traces. This is why we denoted the "hybrid" operators 
(|4,lip by a distinct symbol X (thereby continuing a steady tradition of the field which already 
have Y-systems, T-systems, Q-systems, etc. to occupy all letters of the alphabet). 



5 Functional Relations 

The results of Section [3] imply various functional relations for the Q-operators. To derive them 
we need to use some additional properties of the twist operators (14. 5h which are not immediately 
obvious from their definition (|4.3p . Let Dj and Dj be the operators (|4.5p . corresponding to 
\jj\z) and Lj'(z) from the LHS of (I3.4p , By explicit calculation one can check that the product 
of these operators commuteswith the similarity transformation S defined in (|3,8p , 

[D/X>j,5]=0. (5.1) 

Moreover this product can be rewritten in the form 

V I Vj = V IUJ V G , (5.2) 

where T>j\jj and Vq are the twist operators obtained from (14. 3|) for the operator L/uj(z) and 
the z-independent L-operator G from the RHS of (|3.4p . Again the relation (|5.2p is verified by 
direct calculation, where one needs to take into account the explicit form of (13. lip . (13, 12ft . (13. 13ft 
and flSJID . 

Next, define a scalar factor, cf. (11 .4|) . 

A 7 (*) = A {aiia2) ... iap} (*)= J] 2i Bin ( ** (5.3) 

l<i<j<p ^ ' 

which depends on the set / and the fields $]_, $2j • • • 5 &n- Combining (|3.4p with the definition 
(|4.13p and taking into account (I5.ip and (I5.2P one obtains 

A/X+(z + f , A P1 ) A j X+(z + A - f , A P2 ) = A I{JJ Xj UJ (z, A Pl+P2 ). (5.4) 

6 The Q-operators can be conveniently associated with nodes of a hyercubical Hasse diagram. They are labelled 
by their index sets I = {ai, as, • • • , fflp}, ordered by inclusion. See next section, and Appendix [C] 
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Figure 2: Jacobi determinant formula. 



There are two nontrivial steps in the derivation of the last formula which require explanations. 
First, the simple transfer matrix 

T G = Ttjvm {b g G ® G ® • • • <g> g} = 1 (5.5) 

that arises in the calculations is equal to the identity operator. Second, the scalar factors in (|5.4|) 
arise due to the difference in the definition of the trace over the oscillator algebras (normalized 
trace (|4.10|) ) and over the representation of gl(p) (standard trace). From (|3.15|) it is clear that 
P\Pi oscillator pairs have to be "redistributed" to support the Holstein-Primakoff realization of 
the infinite-dimensional representation ir^ of the algebra gl(p\ + P2)- 
A particular simple case of (|5.4|) arises when P2 = 1 and J = {a p +i}, 

AjX+Cz + I.Ap) Q ap+1 (z + A p+1 -§) = A /Uap+1 X+ Uap+i (z,A p+1 ), (5.6) 

where by the definition (|4.23j) one has Q a (#) = X; a }(z, (0)). Iterating the last formula one 
obtains 

A/ X}(z,A p ) = Q ai (z + \' 1 )Q a2 (z + \' 2 )---Q ap (z + \ , p ), (5.7) 
where the notation here is the same as in (|3.19p and (|3.20p . Next, applying (|4.19p one gets 

A/ Xj(z, A p ) = det|| Q at (z + A$) \\i<i,j< P , (5.8) 

and setting A p = (0) one finally arrives at 

A/ Qj(z) = det||Q ai (^-j + 2±i)||i< <J <p. (5.9) 

Note also that in a particular case / = {1, 2, . . . , n} the formula (|5.8p leads to the determinant 
expression [U El E2 US M\ for transfer matrix KT5h . 

A { i,2,...,„, } T An (z) = det|| Qi(z + A^-) ||i<ij< n , (5.10) 

where A{ j is defined in (I5.3P and A^ in (|3.19p . 

As previously mentioned the 2 n different operators Q/ can be assigned to the nodes of 
a hypercubic Hasse diagram. We will now show that four Q-operators belonging to the same 
quadrilateral as in Fig. [5] satisfy a remarkably simple functional equation, which can be identified 
with the famous Hirota equation from the theory of classical discrete evolution equation. Define 
a matrix 

Mij = Q a .(z - j + S±^) , i,j E {0, . . . ,p + 1} . (5.11) 

where {ao, a%,. . . , a p , a p +i} in an increasing sequence of p + 2 integers which contains the sub- 
sequence / = {ai, . . . ,a p }. Denote a = ao and b = a p +i. Let us now use Jacobi's formula for 
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determinants (see Fig. [2]) for the matrix (|5.1ip . Applying (|5.9p for the subdeterminants one 
obtains the following operatorial functional relation_| 

A{ a ,b} QlLtoU&O) Ql(z) = Q/Uo(^ ~ §) Q/Uft(^ + \) ~ QlUb{z ~ \) QlUa(z + \) , (5.12) 

which was already stated in the introduction as (|1.3p . and where A a ^ is given in (jl.4p . Since all 
Q-operators commute with each other the same relation (|5.12p holds also for the corresponding 
eigenvalues. We will use it to derive Bethe equations in the next section. 



6 Bethe Equations 

The connection between the Hirota equations (|5.12p and the Bethe Ansatz equations is well 
understood [Qj [T8l [19| [20l 136"! 116] . Consider the equations (I5.12P at the eigenvalue level. The 
reader might find it useful to look at the examples of Hasse diagrams in the introduction and in 
appendix ICl when following through the following derivation. Let us denote the zeroes of Qiua(z) 
by 4 Ua - Taking z + \ = z^ a and z - \ = z 1 ^ 11 , equation (I5T21) reads 

Q/UaU6(4 Ua " |) Q/(4 Ua - \) ~ Q/Ua(4 Ua " 1) Q/U6(4 Ufe ) , (6-1) 
Q/UaU^ + \) Q/(4 Ua + |) ~ -Q/Ua(4 Ua + 1) Q/Ub(4 U6 ) , (6.2) 

respectively, where a, 6 /. Taking the ratio of these two equations above one obtains 

n _ Q/(4 Ua " I) Q/Ua(4 Ua + 1) Q/UaUfc(4 Ua " I) 



Q/(4 Ua + \) Q/Ua(4 Ua " 1) Q/UaUb(4 Ua + 



(6.3) 



2 ' 



Here I can also be the empty set. In this case we can remove Qz(z) from the equation using 
(I4.2ip . The number of elements in I UaUb cannot exceed n, therefore / contains at most n — 2 
elements. Thus one obtains n— 1 different relations of the type (|6.3p . with various cardinalities of 
the set /. This exactly matches the number of levels of nested Bethe equations for the g[(n)-spin 
chain. 

Let us take any sequence (a\, . . . ,a n ) of elements of the set {1, . . . , n}. We construct a 
sequence of ascending sets = Iq C I\ C ... C I n = {1, . . . , n} such tlialo Ii = U Oj. Then 
for each Jj, i = 1, . . . , n — 1, we can rewrite (|6.3p as 



(6.4) 



We will call an equation with z 1 * roots the i-th level equation. One need to use (|4.2ip for the 
lowest level, i = 1, and (|4.24|) for the highest level equation with i = n — 1. 

It is not difficult to see that the system of functional equations (|6.3p already corresponds to the 
nested system of Bethe equations of the gl(n) compact magnets. To recover their traditional form, 
we merely need to substitute the eigenvalue formula ([4.25P into (|6.4p . The crucial point of our 
approach is that the analytic structure of (I4.25|) is rigorously derived without any assumptions. 



7 More general relations involving operators (|5.8[) can be obtained in the same way by replacing the arguments 
of the Q-operators in (|5.11[) with arbitrary constants Zj, j = 0, . . . ,p + 1. 

8 This construction can easily be obtained by choosing a path on the Hasse diagram leading from to 
{1, . . . , n}. See Appendix [Cl for more details. 
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It immediately follows from the explicit construction of the Q-operators in Section 4. Then one 
can rewrite (|6.4p as 




1 ' 



(6.5) 



for the lowest level, 

-11 A-i , i lL z i i _ z ii_ 1 LL z Ii_ z i i+1 + l ' vo-dj 

fe Z « Z fc + 2 k& Z l Z k 1 k Z k 2 ^ 2 
for i = 2, . . . , n — 2 and for the highest level 

( In-! 1 \ L I n -1 _ I n -2 _ 1 Zn-1 Z„_i 

Z l + 2 = TT z l z k 2 TT z < z k + 1 /g ^ 

/„_! 1 11 /a-n-1 _ In-2 j_ 1 1 1 _/n-l _ Z„_l _ -, V ' J 

Z l 2 / k Z l k 2 fc^i Z « Z fc - 1 

Equations (|6.5p - (|6.7p are exactly the Bethe equations for the compact gl(n) symmetric spin 
chain. It is obvious that there are n! alternative forms of the above Bethe Ansatz equations, 
corresponding to n! permutations of the elements of the set J, which in turn are associated with 
the n! different bottom-to-top paths on the Hasse diagram. 

To conclude our new solution procedure for the g[(n)-spin chain we give the expression for 
the eigenvalues of ([Lip . It only involves the roots z 1 "- 1 on the last-level 



The derivation of the energy (16.8h from the functional relations is given in Appendix [Dj 

7 Conclusions and Outlook 

In this paper we developed a novel, systematic procedure for constructing the Baxter Q-operators 
connected with the g[(n)-spin chain ([Lip with quasi-periodic boundary conditions. For illustra- 
tion purposes we confined ourselves to the case where the quantum space is a multiple tensor 
product of the compact, fundamental representation of su(n). The Q-operators are constructed 
as transfer matrices associated with infinite-dimensional representations of the Yangian Y(g[(n)) 
built from Fock representations of the harmonic oscillator algebra. This involves rather simple, 
but hitherto unknown solutions (|2.20p of the Yang-Baxter equations with first order dependence 
on the spectral parameter. These solutions provide fundamental building blocks for all other 
required solutions via the standard fusion procedure. As a result we derived the full set of 
functional relations, which enabled us to obtain a new algebraic solution of the su(n)-spin chain 
independent of the Bethe Ansatz. 

The construction we have presented in this paper generalizes in a beautiful fashion to com- 
pact spin chains with su(n\m) supersymmetry. In particular, the partonic Lax operators (| 1 . 1 6 [) 
naturally generalize to the supersymmetric case. Details on this will be reported elsewhere |12| . 
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Our methodology also generalizes to more general representations in the quantum space: 
Recall that in this paper we always took it to be a tensor product of L fundamental ra-dimensional 
representations. Particularly relevant is the case of non-compact spin chains. It will be interesting 
to spell out the exact relation between our construction of Q-operators and a rather different 
approach proposed in |37j and a large number of earlier work cited therein. 
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A Fusion of the Canonical L-operators 

The calculations which one has to perform in order to prove formula (I3.4p are tedious but rather 
straightforward. It is enough to calculate the matrix product of the operator-valued matrices 
(13. 6p and (13. Zj) , and subsequently use the similarity transforms (I3.9|) and (I3.10p . keeping in mind 
the commutation relations between the various oscillators. For the convenience of the reader we 
present here the explicit action of these similarity transforms S\ and 1S2 on the oscillators. 

Similarity transform S\ 



1 ab 1 ab ab 




(A.l) 
(A.2) 



All other oscillators are unchanged. 



Similarity transform S 2 




(A.3) 



cglUJ 




(A.4) 



(A.5) 



All other oscillators are unchanged. 
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B Proof of the Factorization Formula 



In this appendix we prove that the ordered product of the n different partonic Lax operators, 
each containing n — 1 pairs of oscillators, indeed can be disentangled by a similarity transform Sc 
into the product of the st(n) invariant Lax operator C + and a compensating matrix Gc- Further 
details were discussed in section [TJ Instead of using the prescription developed in section [3] to 
derive equation f|l . 18|) inductively, we consider the following decomposition of the partonic Lax 
operators: 

A \ 



L a (z a ) = U" 1 



1 



D a+ i j(l 1 



U a +1 , 



V D n , a 1 / 

with z a = z + A' a . 

The matrices U a do not depend on the spectral parameter. They are of the form 

u«=(i+ E (^EE^f 1 - E h l^X 

\ Kb<c<a ) V 6=1 c=a ) \ a<b<c<n / 



(B.l) 



(B.2) 



The entries in the column below the spectral parameter of the middle matrix on the right hand 
side of equation (IB.lj) take the values 



Dab = -b^ 6 - b ab + E b L b cfe + E h ac h cb ; (B.3) 

c=l c=a+l 

the remaining non-diagonal entries vanish. 

From this decomposition it is obvious that the product of the n partonic Lax operators can 
be written as 



LiOi) • • • L n (z n ) = U 1 



-i 



D 



2.1 



+1) 



(B.4) 



\ Dn,l ' ' ' -D n n _i Z n 2 J 
where two matrices Ui and U n +i only contain an upper triangular part: 



A "b f 12 ••• -b f lrt \ (\ bia ••• bi„ \ 1 



Ui 



b 



V 



t 

n— 1 n 
1 / 



u 



n+1 



\ 



'■ b n _i, 
1 



(B.5) 



Up to now, the lower triangular matrix in equation (|B.4p contains all n — 1 oscillators. A short 
calculation shows that half of them can be absorbed by the similarity transform Sc = S n - ■ ■ S\ 
with 



S a = exp 



a-l 



a-l 



E( b L+ E b ^ ] b i 

c=6+l 



ab 



.0=1 



(B.6) 
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Q{1,2} 




Figure 3: Hasse diagram for gt(2). 

This transformation coincides with the one that can be obtained from the method^ established 
in section [3j Furthermore, it leaves Ui invariant and its action on U n +i factorizes as 

S £ U n+1 SZ 1 = \Z[ 1 V n +i. (B.7) 

As a result one obtains that the ordered product of n partonic Lax operators can be written 



as 



Li(zi) • • • L n (z n ) = ScUx 



\ 



D 



2,1 



where 



D ab = -Kb + ^2 b ac b c6- 



n-1 



c=a+l 



We identify the s[(n) invariant Lax operator in its factorized form (see e.g. [30 

D 2 ,i 



and 



G'C = U n+ i 



(B.8) 



(B.9) 



(B.10) 



(B.ll) 



C Hasse Diagrams 



The quadratic functional relations (j5. 12|) possess an interesting graphical interpretation, cf. Fig.[5j 
The full set of functional equations is then nicely depicted by so-called Hasse diagrams, cf. |16j 



9 As we are interested in the factorization formula where the compensating matrix appears on the right hand 
side of the Lax operator it is convenient to start from the partonic Lax operators L n (z„) and multiply the 
remaining ones in the desired ordering from the left. 
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Q{1,2} 



Q{1, 2,3,4} 




Q{1,2,3} Q{1,2,4} Q{1,3,4} Q{2,3,4} 




Q{3,4} 



Figure 4: Hasse diagram for 0t(4). 

and references therein. These diagrams are used to represent partially ordered sets. In our case we 
take a power set of {1, . . . , n} with order given by the inclusion relation, namely A < B <4> A C B. 
Then Q/ will inherit this ordering, giving us a partially ordered set containing all Q-operators. By 
way of example, Hasse diagrams for n = 2, 3, 4 are presented in the figures E1Q] and 01 respectively. 
To read off the functional relations it is enough to take any 4-cycle in these diagrams, using the 
equivalence depicted in figure [5j 



IUaUb 



QlUa 




IUb 



2isin 2* Q IUaUb (z) Qj(z) = Q IUa (z - |) Q/ufe(^ + \) 

~ QlUb(z ~ \) QlUa(z + |) 



Figure 5: Graphical depiction of the functional relations (|5.12p 

Every path in the Hasse diagram which leads from Q to Q{i.... jT1 } defines a system of 
equivalent but distinct nested Bethe equations. To find each such system, it is enough to take 
all Q-operators on a given path and write one relation for any three subsequent functions on the 
path. Such relation can be written for every three subsequent Q's because there always exists 
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a unique 4-cycle containing them. Finally, it is interesting to point out that the Hasse diagram 
corresponding to the Qt(n) algebra forms an n-dimensional ordered hypercube. 



D Energy Formula 

Let us choose a path on the Hasse diagram given by the ordered set {a±, . . . ,a n } = {1, . . . , n}. 
We define I k = I k -i U {a k } with Iq = 0. The formula for the energy in terms of the Bethe roots 
depends only on the last-level Q's and is given by 



1 

i -. 1 ( _,-*Tl — 1 \ 



fe=l 4 - |^ fe 

where mj n _ 1 is the number of roots of the Q/„_ 1 (^) function. To prove (|6.8p it is enough to take 
the functional relation 

Xl k (zQ, *2, ■ ■ ■ , Zk) _ Xj fc (z , gj, . . . , Z fc _i) X Jfc _ 1 (z 2 , ■ ■ ■ , Zfc) Xj fc _ 1 (z ,Z 2 , • • • ,Zk-l) 

X Ik (z 1 ,z 2 ,...,z k ) X Ik (zi,z 2 ,...,z k ) X Ik _ 1 (zi,z 2 ,...,z k -i) X k _ 1 (z 1 ,z 2 ,...,z k -i) 

(D.2) 

where X/ fc is an eigenvalue of X/, defined in (|4.11|) . It holds for k = 1, . . . , n and originates from 
the fact that every Xj. can be written as a determinant of partonic Q-functions as (see (j5.8j) ) 

A Ik X Ik (z 1 , . . . ,z k ) = det Qa(z b ). (D.3) 

a,bel k 

In this case the relation (1D.2[> is just a version of the Pliicker relations and is valid for any 
numbers z a . In (|D.2|) we omitted an explicit dependence on the spectral parameter which can 
be recover for every Xj k just by taking the mean value of all z's present there. 

If we take k = n and put z a = z — a + ^i— then using this relation recursively we get 

+ j) = y Qv±(f + i + Qin-fc-iCg - j + 1) fD4) 

where Tn denotes an eigenvalue of the transfer matrix in the fundamental representation. In 
order to prove (|6.8|) it is enough to use the well known relation 

E = 2L-2^-logT n (z + -) (D.5) 

az n 2=0 

and the fact that Q/ n (z) = z L where L is the length of the spin chain. 
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